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Abstract 



Fractalization of torus and its transition to chaos in a quasi-periodically forced 
logistic map is re-investigated in relation with a strange nonchaotic attractor, with 
the aid of functional equation for the invariant curve. Existence of fractal torus 
in an interval in parameter space is confirmed by the length and the number of 
extrema of the torus attractor, as well as the Fourier mode analysis. Mechanisms 
of the onset of fractal torus and the transition to chaos are studied in connection 
. with the intermittency. 

"5 : 1 Introduction 

Transition from a torus to chaos has been intensively investigated 0. There two types 
of instabilities exist, which cause the collapse of tori: One is the instability along the 
phase direction and the other in the amplitude direction. The former instability has 
already been studied in detail using the circle map while the oscillation and fractalization 
^ ■ of tori have been reported due to the instability along the amplitude direction. One of 
the authors (KK) reported that the oscillation of tori gets stronger with the increase of 
external forcing, until it reaches fractal, and then chaos appears H . Although the fractal 
nature of the torus was confirmed at the onset of chaos, the strange oscillation of torus 
before the onset remained unclear. 

On the other hand the existence of strange nonchaotic attractors(SNA) was shown 
by Grebogi et al||. Here the word "strange" refers to the geometry of the attractors, 
and the word "chaotic" refers to the orbital instability of the dynamics. In a model 
with quasiperiodic forcing similar to that for the fractalization of torus, they showed 
analytically that the attractor has a non-positive Lyapunov exponent but a complicated 
geometry. Indeed the attractor is not smooth on the set of dense points. Although 
the original model by Grebogi et al. |§ excludes the possibility of chaos, the SNA has 
generally been observed in a system where chaos appears with the increase of the forcing 
as was discussed in the fractalization of torus @, ||, [7], ||, |], p~0[1 . 

In the present paper, we study how the transition from smooth to fractal torus occurs. 
We have confirmed that there is a parameter region with a non-zero measure where the 
SNA exists between the smooth torus and chaos. This leads to the following scenario of 
the destruction of tori. 

Torus— ^Fractal Torus(SNA)— >Chaos 



Figure 1: Three types of attractor: The value of e in each figure is (a)e = 0.1, (b)e = 0.156, 
(c)e = 0.18, respectively. We have plotted 30,000 points after discarding initial transients, 
where a = 3.0, u = (y/E — l)/2. 



Figure 2: Lyapunov exponents are plotted versus the parameter e. Computed from the 
average over 10 6 iterations of the map (1), after discarding initial 5,000 points of transient. 

In the following sections, we will first characterize the nature of fractal torus as a SNA. 
Besides the direct simulation of the quasi-periodically forced logistic map, we adopt the 
functional map to obtain the invariant torus. The length of the invariant curve, as well 
as the number of singular points, shows how the fractalization occurs with the parameter 
change. In §3, perturbation expansion of the functional equation is given, where the 
SNA is characterized by the loss of the convergence in the Fourier mode expansion of 
the derivative of the invariant curve. Comparison between smooth and fractal tori is 
given in §4, where some dynamical signatures of the fractal torus is discussed, such as 
the parameter sensitivity of the Lyapunov exponent, and the phase sensitivity. The onset 
of chaos from the SNA is investigated in §5, where the transition is associated with the 
intermittency from the fixed point in the functional map. 

2 Fractalization of Torus as SNA 

Transition from torus to chaos has been studied with a variety of 2-dimensional 
maps0. With the change of the bifurcation parameter, the amplitude of oscillation 
gets stronger, accompanied with phase- lockings. To focus on the amplitude instability 
by excluding the phase locking, we choose the following two dimensional map [[|, with a 
constant rotation of the phase; 

f x n+1 = f(x n ) +eg(9 n ) m 
t n+1 = n + u;mod 1 1 J 

Here x n represents the amplitude, while 9 n corresponds to the phase of oscillation. To 
exclude the phase lockings no coupling from x to 9 is included, where uo is an irrational 
number representing the rotation number. 

In this paper we set f(x) = ax (I — x),g(9) = sin27r6', and to = (y/E — l)/2, as the 
simplest nonlinear map with an external driving force. For most simulations, we fix 
a = 3.0, and vary e to see the change of the attractor. At e = 0, the map, of course, is 
just a one-dimensional map of x n with a fix point x* = 1 — 1/ a. Thus the attractor of the 
two-dimensional map is just a with a straight line torus. As e is increased, oscillations 
of the torus start to appear, which become fractal with the dimension above one ( as 
is confirmed later by direct measurement of the length of the torus). As e is increased 
further, another transition from fractal torus to chaos appears, by which the stability 
along the x-direction is lost. 

Three examples of the patterns of attractors are given in Fig.|l] corresponding to the 
three types of attractors, while the change of the Lyapunov exponent is plotted in Fig.0. 
Transition from torus to chaos occurs at e ~ 0.1573. As will be confirmed later, the 
torus loses its smoothness around e ~ 0.1553. Hence the SNA exists in an interval of the 



Figure 3: Examples of log-log plot of length versus mesh (a = 3.0, F k = 317811), obtained 
from the attractor of the functional map. 

parameter, where the Lyapunov exponent has a sharp sensitivity to e, as in the chaos 
region. 

To confirm that the attractor has a non-integer dimension, we study the equation 
for the invariant torus, following the argument of ||. In our model (|I|) the attractor is 
expressed as a single- valued function of 9, as x = X{9) ( < 9 < 1 ). If this function 
represents an invariant curve of the map, it must satisfy the following functional equation. 



X(fl + cj;modl) = f(X(0)) + eg(6) 

= aX(6)(l-X(6))+esm27i6 (2) 

This equation is postulated by the constraint that the point (x n+ i,9 n+ i) should be 
also on the curve, and is obtained by substituting x n by X(9) and x n+ \ by X{6 + uj) in 
the map ([!]). From this functional equation we can introduce the following "functional 
map" , which maps a function to another function: 



X n+1 (# + ^;modl) = f(X n {9)) + eg(9) 

= aX n (6)(l-X n (6))+esm2Tr6. (3) 

The attractor of the original map (|I|) is obtained as a fixed point in the functional 
space, for the iteration of the functional map (|3]). This discussion stands only if the 
Lyapunov exponent is negative and the attractor has a single value x for each 9. 

Since the functional map (|^) is an infinite dimensional dynamical system, one cannot 
compute it directly. We have numerically computed it by approximating to = (y5 — l)/2 
by u k = F k _i/F kl where {-FA;}fc=o,i,... are Fibonacci series. {oj k — > u as k — > oo ). This 
approximation transforms the functional map to an F k - dimensional map which maps F k 
lattice points on the 9 coordinate onto themselves. We have computed the attractor of 
this Ffe-dimensional map, to obtain the approximate solution of the functional equation(|2]) 
as a piecewise-linear function. 

For the parameter regime corresponding to the torus attractor, the convergence is 
rather fast. The convergence time of the functional map to a fixed point gets longer as e 
approaches the onset of chaos, while, for e > 0.1573, corresponding to the chaos region, 
the functional map does not converge to a fixed point. 

The functional equation enables us to compute the length Lj of X(9) as follows: 



F k -l 



If the dimension of the attractor is 1 + a, Lj must scale as Lj oc j~ a , where < a < 1. 
The slope of this plot a gives a kind of fractal dimension. Fig.|3] is an example of the 
log-log plot of length Lj versus the mesh width j. Here we have adopted F k = 317811 for 
the approximation of uj. Hence j = 1 corresponds to the mesh size of 1/317811 ~ 3 x 10 -6 . 
Fig.|| is the plot of slope a versus external force e. 

The slope starts to be nonzero at the transition point at about e ~ 0.1548. The slope 
jumps to a ~ 0.62, and stays around the value with the increase of e, up to the vicinity 
of the onset of chaos. 



Figure 4: Change of the slope a versus external forcing s (a — 3.0, Fk = 317811). The 
slope is estimated from the length at the 5 smallest mesh scales. The error bar is not 
explicitly given, but it is about .05 ~ .1. 



Figure 5: Length of attractors with the mesh scale j. F k = 9227465, and a = 3.0. 

In the region between e ~ .1548 and .1553, the exponent is roughly 0.3 up to fine mesh 
scales (e.g., 1/317811), but there appears a saturation at a finer mesh. In Fig.|5| we have 
plotted the length versus mesh, taking F k = 9227465. There is a crossover to a smooth 
behavior at the scale about 10/9227465 for e = 0.1548. Here we call this parameter regime 
as "pre-fractal-torus" , since this region is distinguished from the smooth and fractal torus 
regimes. Indeed the crossover to a smooth curve is seen only at a much finer scale than 
for a smooth torus regime, while this crossover scale increases as e approaches the onset 
parameter for fractal torus. The slope at the scaling regime (i.e., scales with a larger 
mesh beyond the crossover) has a jump at the onset of fractal torus: The exponent a in 
the scaling regime jumps from .3 to .62, when the parameter e crosses 0.1553, the border 
between pre-fractal and fractal torus regimes. 

For e > 0.1573, X n (9) does not converge to a fixed curve. Indeed the Lyapunov 
exponent of the two-dimensional map is positive in this regime, and the attractor here 
is chaotic. Summing up the results of the fractal exponent of Fig.£|, and the Lyapunov 
exponent of Fig.^], one can conclude that the fractal torus exists at least in the region 
0.1553 < e < 0.1573 l . 

We have also measured the distribution of local dimensions, by varying e, which is 
shown in Fig.^. The distribution is single-humped, with slight asymmetry. At the "pre- 
fractal-torus" , the distribution has a peak around zero. 

We have also calculated the number of extremal points of the function X(9). Here 
the number Nj of extremal points for a given mesh is defined as follows. First count the 
number of the mesh point % such that 

x (^)- x (±)\t x (i±li)- x (i±i)\ <0 (5) 



Fk ) \FkJ J I \ Fk J \ Fk 

is satisfied: Then Nj is obtained by dividing this number by j. Fig.[?| gives some examples 
of the log-log plot of Nj versus j. 

The number Nj increases as the mesh j is reduced, which means that the finer mesh 
one uses, the more extrema appear. This is consistent with the fractal nature of torus. 
However, we note that Nj scales not as Nj ~ j -1 but as Nj ~ j _/3 , where < (3 < 1. 
This means that the fraction of singular points decreases with the mesh. In other words, 
the extremal points lie on a Cantor set on the 9 axis. 



1 By numerical means, one cannot completely exclude the possibility that the the non-smoothness of 
the torus disappears as the mesh scale goes down to a much finer scale, and the fractal torus exists 
only at the onset of chaos. Indeed this suspect has lead to the report of the fractal torus only at a 
single parameter point, in We assume that the scaling with the exponent a lasts to an infinitesimal 
scale, since the increase of the crossover size near 0.1553 is so strong that it is distinguished from the 
pre-fractal-torus region. 



Figure 6: Distribution of the local dimension obtained through the measurement of the 
local length. We have measured the length of X{9) at each lattice point using 6000 lattice 
points around it. By computing the local dimension over all lattice points, and sampling 
the number with the bin size 0.01, the distribution is obtained, (a = 3.0) 



Figure 7: The number of extremal points versus mesh size. See text for the method of 
computation, (a = 3.0, F k = 317811). 

Figure 8: (a),(b)Change of the Lyapunov exponent with the parameter e(a = 2.8, u = 
— l)/2). (b) is the blowup of (a) for 0.204 < e < 0.208. The exponent is computed 
from the average of 10 6 steps after discarding initial 5000 steps. 

Dependence of (3 on the parameter e is plotted in Fig.f|. The exponent stays around 
0.75 while the fractal dimension stays around 0.6. Near the onset of chaos, e = 0.1573, 
the exponent (3 increases, till it takes 1.0 at the onset. In other words, the ruggedness 
of the torus fills the space down to any small scale when the torus loses stability and is 
replaced by chaos. 

The linear increase of the extremum point agrees with the picture by a random curve. 
Consider, for example, a curve C(i) (with i lattice point) generated by C(i+1) = 5xC(i) + 
7) with rj a random number distributed over [-1,1] and 5 < 1. The number of extremum 
points increases linearly with the number of mesh points. The increase is given by 2iV/3 
for 5 = 0, while the increase is slower with a smaller proportion coefficient for 5 > 0, till it 
is given by N/2 for 5 = 1. Since the increase is given by 0.58iV at the onset of chaos (and 
also in the chaos region) in our simulation, the critical torus can be approximated by a 
correlated random curve. This behavior and the coefficient characterize how the fractal 
torus collapses and is replaced by chaos in the two dimensional map eq. ([]]). 

To close the present section, it is interesting to note the universality of our results. In 
the discovery of a fractal torus by [Q, the parameter a = — 1 was adopted, where successive 
transition from torus to fractal torus, and then to chaos was found with the increase of 
e. We have re-examined the simulation using a larger mesh sizes. Here again, a fractal 
torus exists for a finite interval of parameters, where the fractal dimension increases from 
0.6 to about 0.88 with e. The number of extremal points again increases with a fractional 
power, while it increases linearly at the onset of chaos. 

Another interesting set of examples is given for a = 2.8. The Lyapunov exponent is 
plotted in Fig.|, while the fractal dimension and the number of extremal points obtained 
from the functional equation (with the mesh 317811) are given in Fig.|S|. 

From these figures, one can see that there appears windows of fractal torus (FT) 
and smooth torus(ST), beyond the onset of chaos. Indeed there are transition sequences 
among ST, FT, chaos(C), as is seen in Fig.|8| and in Fig.|9]. At e ~ .204 FT appears, which 
is replaced by chaos at e ~ .206. For e ~ .208, the "inverse" bifurcation from chaos to 
FT, and then to ST (at e ~ .211) proceeds. Further bifurcations to FT and chaos, and 
back to ST are seen for larger e. Again at the boundary between chaos and fractal torus, 
the increase of the extremum points is almost linear. 

3 Fractal Torus viewed from Functional Map 

In this section we apply the Fourier mode analysis to the functional equation(fJ). First 

Figure 9: Fractal dimension a (obtained by measuring the length) and the exponent (3 
on the number of extremal points are plotted versus e. Computed from the attractor of 
the functional map. (a = 2.8, F k = 317811). 



we consider the Fourier expansion of X(8) 



X(6) = £ X{k)e 2 ™ ke (6) 

k=— oo 

= \ X(9)e- 2mke d6 (7) 



o 



and substitute it into @. Then we have 

X(A;)e 2 ^ = aX(fc) - a £ X(k')X(k - k') - —(S k ,i - (8) 

k'=—oo ~ 

Next we expand each Fourier mode with respect to the powers of e. 

oo 

X{k) = Y,e n X n {k) (9) 

n=0 

Substituting this into (§) and comparing the terms of same power of e in both sides of 
equation, we get 



oo n 



X n {k)e 2 ^ = aX n {k) - a £ X m (k')X n „ m (k - k') - ~5 n ,i(4,i - 4 ,-i) (10) 

k'=— oo m=0 

for each n. It is straightforward to show that has only terms e n such that n > \k\. 
This means that the term e n only has a Fourier mode with the wave number no more 
than n. The first three terms of X n (k), X(k) and X(9) are given as follows: 

Mk) = (l-~)*k,o (11) 

= 2(e^+a-2) (12) 



^/.x = -a{2X 1 (-l)X 1 (l) ( 5 M + (^(-l)) 2 ^ + (^i(l)) 2 4, 2 } ( , 



1(0) = X (0)+£ 2 X 2 (0)+£ 4 X 4 (0) + --- (14) 
X(±l) = £ X 1 (±1)+£ 3 X 3 (±1)+£ 5 X 5 (±1) + --- (15) 
X(±2) = e 2 X 2 (±2)+£ 4 X 4 (±2) + --- (16) 



X(0) = X o (0)+e{Xi(l)e^ + Xi(-l)e-^} 

+ £ 2 {x 2 (0) + X 2 (2)e 4me + X 2 (-2)e" 4 ^} 

+e 3 {X 3 (l)e 2 ^ + X 3 (-l)e- 2 ^ + X 3 (3)e 6 ™ e + X^(-?,)e~ &m6 ) + • • -(17) 

When e is gradually increased from zero, the terms of higher frequency modes get 
larger in order. This explains the amplification of torus oscillation for larger e, as well as 
the slower decay of the Fourier coefficient with the wavenumber. 

We have numerically calculated the power spectrum P(k) = \X(k)\ 2 for the three types 
of attractors, by using the largest 2^ points possible out of F k points of the attractor of 
functional map ([!]). The power spectra P{k) are plotted in Fig.|l0|, for a = 3.0,e = 
0.15,0.156,0.16, where F k = 317811, 2 N = 2 18 = 262144. 



Figure 10: Power spectrum P(k) = \X(k)\ 2 for a = 3.0 and (a)e = 0.15, (b)e = 0.156, 
(c)e = 0.16 obtained by averaging 100 different phase shifts 9 (shift of lattice points 
within the range of 1/F k ). (F k = 317811, 2 N = 2 18 = 262144) 

When the invariant curve is fractal, the first derivative of the Fourier series^) is 
expected to lose its convergence. Indeed, the spectra P(k) decrease with k slowly for 
fractal torus. The maximum of the envelope of P(k) decays slower than k~ 2 , which 
means that X(k) ~ k~ a , where < a < 1. The non-smoothness is seen in the power 
spectra. For a smooth torus the spectra decay faster than or equal to k~ 2 . 

It is interesting to note the analogy with the cascade process in the turbulence. The 
equation (||) has a similar form with the Fourier expansion of the Navier-Stokes equation, 
with regards to the formation of cascade process by the term X(k)X(k — k'). In the case 
of fluid turbulence, this term brings about the Kolmogorov's energy cascade. In contrast 
with the time-dependence in the turbulence, the Fourier modes of our functional equation 
fall into time-independent values at the fractal torus regime. In the fluid turbulence, the 
intermittency leads to sporadicity in the vortex cascade, which brings about the deviation 
from Kolmogorov's 5/3 law. In our problem this corresponds to the sporadicity as to the 
extremal points. The degree of the sporadicity changes with the parameter e, as in the 
intermittency effect in the turbulence, while the sporadicity is lost at the the onset of 
chaos in the cascade of the functional equation. It will be interesting to search for a 
quantity in our problem corresponding to the energy cascade in turbulence. 



4 Fractal Torus versus Smooth Torus 



Let us reconsider the difference between smooth and fractal tori from the dynamical 
systems viewpoint. Of course, the largest difference between the two types of tori is the 
length. The length of fractal torus is infinite, in contrast with the smooth torus. This 
also leads to a difference in the orbital instability for the dynamics of the two types of 
tori. Smooth and fractal tori are both stable against the perturbation in the x-direction. 
Although no exponential divergence of orbits exists for both, they are different with each 
other as to the phase sensitivity. Pikovsky and Feudel j|] have shown that two points 
on the SNA with close ^-values separate from each other, by introducing the following 
phase sensitivity exponent: For this, note that the absolute value of the first derivative of 
the orbit |%H fluctuates with time and sometimes has a large burst. An arbitrarily large 
burst can appear when the map is iterated for infinite time steps. By differentiating eq. 
(HI) with respect to 9, one obtains 

^ = f'M^ + e 9 >(e) (18) 

Using this equation one can compute the evolution of ^ starting from some initial point 
(xq, 9q) and %p = 0. The phase sensitivity by Pikovsky and Feudel, then, is defined as 



r^y = min max 

x fio 0<n<N 



8x r 



89 



This quantity must grow infinitely for SNA with some power (and exponentially for a 
chaotic attractor) Q . In our example, this quantity increases with the power 2.5~3 for 
a fractal region, as is shown in Fig.|TT|(a). At the "pre- fractal-torus" regime the quantity 



Figure 11: (a)Tjv and (b) T^/ obtained for a = 3.0, uo = (V5 — l)/2. We have chosen 
100 initial points (xq,9o) randomly and iterated eqs. (18) and (19) for each initial point, 
starting from ^§§- = 0. T N and T$ are obtained as the minimum of these 100 trajectories. 



increases up to 10 6 steps (which corresponds to the mesh scale of 10~ 6 ), and then saturates, 
as is consistent with the results of the length in §2 obtained from the functional equation. 

This T N corresponds to the largest derivative of the attractor X{6) of functional 
map(B[) for the mesh size about Nx (number of samples for min ). The larger the iteration 

20,00 

number N is, the larger the length of the attractor at the mesh is. 

The second difference between smooth and fractal tori lies in the parameter dependence 
of our invariant curve. For a fractal torus, there is sensitivity of the curve to parameter 
e at least for some value of 9. To confirm the sensitivity, we introduce the parameter 
sensitivity, in the same way as Tjv[Q: By differentiating the first equation of (|1|) with 



respect to e, we obtain 



Then we define 



^ = f'M^ + g(9). (19) 



1 at — mm max 

20,6*0 0<n<N 



dx r 



de 

In Fig.|TTl(b) we have plotted this parameter sensitivity. Again this quantity grows with 
some power for the fractal torus. In our model this power is identical with that for T^, 
i.e., the power 2.5~3 (see Fig. |ll|(b)). This is expected since eq . fll8|) for has the same 
form as (eq.([0|)) except the last term that is smooth. The power law increase of T$ 
implies that the fractal torus always has this kind of parameter sensitivity. 

Such parameter dependence is also reflected in the Lyapunov exponent of each attrac- 
tor. The sensitivity of the Lyapunov exponent to the parameter e is relatively sharp in 
the region of fractal torus, while the parameter dependence is smooth in the region of 
normal torus, as we have mentioned earlier. 

The third difference lies in the loss of convergence in the derivative of X{9) for the 
fractal torus. Indeed, the Fourier expansion of X(9) decays faster than or equal to 1/k 
for a smooth torus, while for the fractal one, it decays slower than 1/k, as is discussed in 
the previous section in terms of P(k). 

This loss of convergence corresponds to the power spectrum of the time series of x n , 
since the relation 9 = um + const, (modi) leads to a direct correspondence between 9 
and time. It is shown in |J that the power spectrum of time series is singular for SNA. 
There the number of singular points whose power is larger than a threshold decreases with 
some power with the threshold, while it decays exponentially for a smooth torus. This 
characteristic is another manifestation of the loss of convergence in the Fourier modes of 
the derivative. 

The question remains what causes the transition from smooth to fractal torus. One 
possible route to SNA is crisis between stable and unstable invariant curves 0]. In our 
model, however, the unstable torus (continued from the unstable fixed point of the logistic 
map) lies far away from the stable one and has nothing to do with the fractalization. Hence 
there must be another possible mechanism for the emergence of SNA. In the functional 
equation, this can be seen in the loss of convergence in the Fourier series. How is the 
loss of smoothness expressed in terms of dynamical systems for the original 2-dimensional 
map? 

To study this problem, we have introduced the following i^-map. If we approximate 
uj by Fk-i/Fk and start iterations from any point (xq,9q), the orbit will be periodic in 



Figure 12: Graph of F^-map for F& = 55, 6 = 



Figure 13: Bifurcation diagram obtained from F^-map for F^ = 987. A sequence of x n 
generated by the F^-map is plotted over 100 steps after 2900 points are discarded. The 
time series is plotted for each e-value. Initial condition for the next e value is chosen 
from the orbit x n for the previous e, while e is incremented by 0.0001 between 0.154 and 
0.1585. 

9 with period F^. (0F k = #o + F^uj = 8o + Fk-i = ^(modl)) Thus a composition one- 
dimensional map is constructed as the F^-times iterations of the original 2-dimensional 
map. Let us call this composition map Ft-map. F^-map is a function of x if we fix the 
initial 9q. The shape of graph of F^-map depends on initial 9 , e and k. Fig.[T2| gives a 
sequence of examples of the F^-map. The map converges to the functional map as k goes 
to infinity. We can infer the behavior of the functional map from the F^-map. 

In Fig.[L3]we have plotted the bifurcation diagram of the F^-map, by taking Ft = 987. ( 
This rather small value of Ft is chosen so that the characteristic feature of one-dimensional 
map is visible). Since the onset of fractal torus and chaos can slightly differ between 
lu = Ffc_i/Ffe and uo = (y/5 — l)/2, we have to take into account of possible shift of the 
bifurcation parameter, to compare the result of Ft-map with that of the functional map. 



Let us look at the SNA region (0.1542 < e < 0.1573) in Fig.|13|. First, we can see 
sensitive dependence on e, as we have discussed earlier. Second, we can see some points 
where the map does not approach a fixed point, and instead, shows chaotic behavior. By 
the definition of SNA, the F^-map must have a fixed point when k goes to infinity. When 
uj is approximated by a rational ( say by F^-i/Fk with Ft = 987), this does not have to 
be true. It seems that there always exists at least one chaotic point for some initial 9q at 
every value of e in the SNA region. From this result we propose the following picture of 
SNA viewed from the F^-map. In the SNA region, the F^-map has dense e values leading 
to chaotic behavior for some #o> whose measure goes to zero with the increase of the mesh 
size Ffc. This kind of "partially chaotic" behavior may be related with the existence of 
positive local Lyapunov exponent discussed by Pikovsky and Feudel[|[]. 

We also note that F^-map has strong sensitivity to #o for SNA. Change of Fk(x;9o) 
with #o gets faster as k goes to infinity. This feature of F^-map is consistent with the 
fractal nature of torus. 



5 From Fractal Torus to Chaos 



What characterizes the transition from fractal torus to chaos? First, it is clear that 
the functional map(^) loses its convergence. This is interpreted as the loss of stability 
of the fixed point in the functional space. The functional map shows chaotic dynamics 
( possibly having a high- dimensional attractor) there, although the snapshot pattern of 
X n {9) remains to be fractal with the power 0.88 ~ 1. Beyond the onset of chaos, the 
chaotic fluctuation of \X n+ i{6) — X n {6)\ 2 increases with e. It is interesting to study 
the transition from a fixed point to chaos in the functional space as a high-dimensional 
dynamical system. As for the bifurcation to chaos, this transition might be referred as 
intermittent chaos in the functional space. 

To verify this idea we have again considered the F& map, the composite one- dimensional 
map, given in Fig. 12. From the graph a kind of intermittent transition is seen at about 
e ~ 0.1561. Indeed there is a tangent bifurcation from the fixed point, by which the 



attractor of the Fk- map changes to a chaotic one. This dependence on e does not change 
qualitatively with k, if k is large enough. Then the intermittent character of transition 
from fractal torus to chaos in the functional map is inferred by taking the limit k — > oo. 

As k is increased, the sensitivity on e will be sharper, not only for a fractal torus but 
also for a chaotic attractor. The sensitivity of the shape of graph leads to sensitivity of 
the Lyapunov exponent. This gives an explanation to the sensitive dependence of the 
Lyapunov exponent on e, previously mentioned. 



6 Summary and Discussion 



We have verified the existence of SNA in our quasi-periodically forced logistic map by 
measuring the Lyapunov exponent, the length and the number of extrema of attractor 
and power spectra. We have found that SNA exists in a finite interval in the parameter 
space, and concluded that the SNA is no other than the fractalization of torus in 0. It 
should be noted that Anishchenko et al. || have also recently discussed the fractalization 
of torus to chaos as SNA, by using a forced circle map. 

To understand the mechanism of fractalization and transition to chaos, we have in- 
troduced a functional map (functional equation ). The attractor of the 2-dimensional map 
is represented as a fixed point of the functional map in the functional space. In this 
view point, fractalization of torus is expressed as the change of the fixed point solution of 
the functional map, from smooth to non-smooth one. The transition to chaos might be 
regarded as an intermittent chaos in the functional space. 

We have also introduced i^-map, a composition of the map over steps. This F^- 
map is a kind of cross section of the functional map. With this map, the origin of the 
fractalization of torus is related with the sensitive dependence on the initial phase 9q, 
while the transition to chaos is associated with the intermittency from the fixed point in 
this Ffc-map. 

There has been general interest in search for a novel type of dynamics in a forced 



system [IT]], as has been pioneered by Moser|12[]. There the instability arising from the 



off-diagonal element |13|| in the Jacobi matrix may lead to a novel dynamical behavior. 



Such instability is also seen in the convective instability in open flow problem ||14j|, where 
the instability due to the off-diagonal element leads to a rich variety of dynamics ||15|| . The 



strange non-chaotic attractor in our model is also due to the instability by the off-diagonal 
element, as is characterized by the phase sensitivity of the amplitude. It is interesting 
to explore dynamical systems with such type of instability in general, where the present 
functional map method may be useful. 

Study of the functional equation itself is an interesting topic. In particular, the at- 
tractor of the functional map is high-dimensional chaos, when the fractal torus is replaced 
by chaos. The one-dimensional string X(9) is spatially fractal and temporally chaotic, 
as in the problem of developed turbulence. In this sense, the model may provide a novel 
class of spatiotemporal chaos. Indeed, the functional equation can be viewed as a coupled 
map with a long-range coupling, by which each lattice point is coupled with the point at 
a distance of -F^-i sites. 

It may be interesting to note an analogy between our functional equation and the so 
called Weierstrass function. Weierstrass function is defined as 

oo 

W(9) = J2 an cos b n 7i9 (a6>l,0<a<l) (20) 

n=0 



Topological dimension of Weierstrass function is one but its fractal dimension is believed 



to be more than one. The Fourier coefficient of W(9) decreases as k lo ^ b with the 



wavenumber k. As Yamaguchi and Hata[16| have demonstrated that the fractal nature 



of Weierstrass (and Takagi) functions are derived from the functional equation. Unfor- 
tunately their analysis is not applicable to ours, since their functional equation is linear 
in contrast with ours. Still some analytical studies on the functional equation as well 
as the renormalization group analysis |17]] for the perturbation series (|T0|) should be of 
importance in future. 
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